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Abstract 



Integrating out a heavy field gives rise to effective Lagrangian containing higher dimensional 

operators. In the context of Type-I seesaw mechanism, integrating out the heavy right handed 
MH 1 neutrino field leads to unique dimension five operator which gives the tree level neutrino mass 

term. Apart from these there are dimension six operators that can have important implications. 
QJ ■ A linear combination of two such operators gives rise to the non-unitarity in the lepton mixing 

matrix, E/pmns- In this paper, we discuss the origin of non-unitarity at the high scale and its 

evolution through renormalization group running. 
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1 Introduction 

Standard Model (SM) of particle physics does not admit mass terms for the neutrinos. This was 
motivated by the absence of observational evidence of right handed neutrinos. However, in the 
last decades, neutrino interferometry experiments have observed oscillation of neutrinos. This is 
ascribed to neutrino mass and their mixing among different generations. Neutrino oscillation data 
gives neutrino mass squared differences, Am|i and |Am 2 | , to the order of O (l0~ 5 ) and (lCP 3 ) 
eV 2 respectively pQ. Also from cosmology, we have an upper bound on the sum of neutrino masses 
[2] which is Yl m i — 0-5 eV. From the results of these observations it is clear that the neutrino 
masses are orders of magnitude smaller than other fermions i.e. charged leptons and quarks. A 
standard procedure to generate such a small neutrino mass naturally is to add heavy fields to the 
SM particle spectra. Below the mass scale of the heavy field, it gets integrated out from the theory 
and gives rise to the effective Lagrangian containing higher dimensional operators. These operators 
are non-renormalizable, and SM gauge group invariant. The lowest dimensional of such operators 
is the dimension-5 operator [3]. 

A ~ of (W Li eA (cf> T e T l Lj ) (1.1) 



where II and eft are respectively the SM lepton and Higgs doublets. As we can see from the above 
expression, c^> has the dimension of M _1 , where M is the mass of the right handed heavy particle. 
This operator violates lepton number by two units and leads to Majorana masses for neutrinos. 
After Spontaneous Symmetry Breaking (SSB), the cfi field acquires Vacuum Expectation Value 
(VEV) v and consequently neutrino mass becomes; m v ~ 2 C v 2 - With v ~ 246 GeV, to produce 
neutrino mass of ~ 0.5 eV, M has to be of the order of O(10 14 GeV). This way of generating small 
neutrino masses, by integrating out the heavy field is popularly called seesaw mechanism (Type-I) 

The next set of higher dimensional operators that arise in the theory are the dimension six 
operators. List of different dimension six operators can be found in [5]. A linear combination of 
two such operators is of special relevance in connection with the neutrino mass and mixings [6j. 
This operator is expressed as 

,(6) 



C e ~ c% (l Li <t>) i$ (fl Lj ) (1-2) 

Note that this operator is of Dirac type in contrast with the Majorana type operator (jl.ip . After 
symmetry breaking, the right hand side of Eq. (jl.2p becomes 

4j ] -^VLii0VLj (1-3) 

which makes the coefficient of neutrino kinetic term deviate from unity by an amount of ^ c v 2 ■ 
Now, in order to bring the kinetic term in its usual canonical form, the field vl is rescaled with 

1 /9 

a factor of (l + ^c v 2 ) ■ This also modifies the charged current interaction. Consequently the 
lepton mixing matrix, £/pmns> becomes non-unitary by a factor of (l — \c^>v 2 ) [7]. Because of 
this, the operator in Eq. (j 1 . 2 [) is called the non-unitary operator and c' 6 ' is called the non-unitary 
parameter. The dimension of c^> is M~ 2 . With M is of the order of 0(1O 14 GeV), non-unitarity 
is highly suppressed. 

Signature of such heavy neutrinos cannot be tested at colliders as their mass scale is way 
beyond the reach of these experiments. Therefore, from the point of view of testability at collider, 
one needs to bring down the seesaw scale to TeV. One way to circumvent this problem is provided 
by inverse seesaw [8]. In these scenarios {i.e. TeV scale seesaw [U [TO]), mass scale of the 
heavy field is of the order of O(10 3 GeV) and hence non-unitarity can be appreciable i.e. c^'v 2 ~ 
0(1O~ 2 ). Consequently the deviation of C/pmns from unitarity can be significant and can have 
testable consequences [TT]. 



As we know that these type of operators arise at comparatively higher scale whereas observations 
are made at a much lower scale [12] , the effect of Renormalization Group Evolution (RGE) needs 
to be incorporated. RGE for the dimension five operator, in the context of seesaw mechanism, is 
well known in the literature |13|, 1X4] . RGE of the non-unitary operator is, however, not as widely 
studied. In [6] only part of the /3-function of the non unitary parameter is computed. Here we 
shall present a detailed analysis of the origin of the non-unitary operator in the context of Type-I 
seesaw and one loop correction to that. Since the non-unitary operator becomes significant in TeV 
scale models, the RG running of this operator from the TeV scale to the electro- weak symmetry 
breaking scale is not expected to be as large as operators which are significant at very high energy 
(10 12 — 10 14 GeV). But as neutrino physics has entered the precision era, so strictly speaking, RGE 
of the non-unitary operator should be explored. Also for some cases, e.g., where symmetry enforces 
C/e3 to be zero, a small RG running could evade the symmetry and may give a non-zero value of 
U e 3 as indicated by present data pp. 

The plan of the paper is as follows; in section 2 we shall present the Lagrangian, origin of 
non-unitarity at high scale and its renormalization procedure. In section 3 we give the one loop 
correction to the non-unitary operator, lepton doublet and Higgs doublet. Section 4 contains 
/3-function of the non-unitary operator. We summarize in section 5. Appendix A contains the 
Feynman rules and appendix B contains the divergent part of different Feynman diagrams. 

2 Origin of non-unitarity and Renormalization procedure 

We add right handed gauge singlet heavy field to the SM particle spectra with mass is Mj. The 
extended part of the Lagrangian is 



C 



Yuk 



^Mass 



±M ij N?P R N j + )i.c. 



(2.1) 



where <j) denotes SU(2) higgs doublet with <$ = ecf)*; e is the 2x2 anti-symmetric tensor. The 
lepton doublet of generation i is denoted by k, and Nj denote the heavy neutrino fields. The 
Yukawa coupling Y v is complex-valued % x j matrix in general. After SSB one gets the Dirac mass 
matrix for the neutrino as mo ~ vY v . The Major ana mass matrix M is j x j complex symmetric 
matrix. 

2.1 Origin of non-unitarity 



With the above Lagrangian (12. ip one can construct four possible dimension-6 Dirac type operators 
in the context of l<p — > 1(f) and 1(f)* — > 1(f)*, two in each cases. Here we present one from each of them. 
The other two will be the same with fermion lines reversed. 
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Figure 1: Feynman diagram, l<f> — > l<j> mediated by heavy right handed neutrino (on the left hand side). On 
the right hand side, low momentum transfer approximation, leads to c-®' vertex. 



Tree level matching of both sides, in the limit s <C M|, gives 



(y u m- 1 (m- 1 )* : yj 



j i 



e-y^lft+fy 



Pl 



l _ (6) 



+ 



Pl 



(2.2) 



We identify c (6) = 2Y i ,M~ 1 (M" 1 )* Fj. Note that c (6) is Hermitian. Here s = (pi + p<j>) 2 . A second 
one is as follows 
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Figure 2: Feynman diagram, 1(f)* — » /</>* mediated by heavy right handed neutrino (on the left hand side). 
On the right hand side, low momentum transfer approximation, leads to c' 6 - 1 vertex. 



Tree level matching of both sides, in the limit t <C M?, gives 
(y u M- 1 (M- 1 )* yj) e 7/3 e Q(S (j*, -^) Pl 



« (6) 



a5(^-/^)ft (2-3) 



Here £ = (p/ — p.) . Note that the external scalar lines as well as the momentum dependence in 
s-channel and £-channel are different. For the dimension-5 operator both s-channel and £-channel 
contribute to the same effective vertex. But here they contribute to different vertices. The effective 
dimension-6 operator (at low energy) in the Lagrangian is given by 



£ (6) = ~4? (W) {4??Iia) 



(2.4) 



2.2 Renormalization procedure 

Here we consider the effect of incorporating one loop corrections to the non-unitary operator (j2.4|) 
i.e. renormalization of this operator. This requires wave function renormalization of II and (ft 
doublet as well as the c^' vertex renormalization. The wave function renormalization is defined as 



{hi) 



B 
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z 2 



I, 






Zh 



(2.5) 



The subscript B stands for bare a quantity. Fields on the right hand side are the renormalized 
quantities. We shall follow the dimensional regularization procedure. Z's are the renormalization 
constants which contain the divergence arising from quantum corrections. Among different possi- 
bilities, we shall follow minimal subtraction (MS) scheme. In this scheme the divergent part of a 
loop diagram is absorbed in the counter term, defined as 






(2.6) 



where e is 4 — d and 5Z j are coefficients of - in the divergent part of a loop diagram. In the same 
spirit, the bare Lagrangian can be written as the sum of renormalized Lagrangian and the counter 
term 



£(6) _ £(6) + ^(6) 



with S£^> , given by 



J^ 6 ) = \bj$^ep)0[f?\j>l 



Ac*- 6 ' is the counterterm, related to the bare and renormalized coupling as below 



.(6) 



VV c(6) + Ac(6) ^i? z 



(2.7) 



(2. 



(2.J 



where \x is the renormalization scale. 

3 One loop correction 

In this section we shall present the one loop correction to the lepton and Higgs self energy as well 
as the c' 6 ) vertex. Here we shall give only the one loop diagrams and the corresponding counter 
term. Expressions for individual diagram can be found in appendix B. In the MS scheme, the 
counterterm is determined by compensating the divergence arising from all possible one loop 1 
particle irreducible diagrams. Pictorially it is shown below for c^> 



+ ^One-loop 1PI diagrams 




j //Pl P$ 





UV finite 



Figure 3: One loop renormalized c^' vertex (on the left hand side). First diagram on the right hand side 
is the tree level c* 6 ' vertex. Diagram on the extreme right is the counterterm vertex which eventually cancels 
the divergence. 

the above figure can be expressed in terms of words as 

Divergent part of all diagrams + the counterterm = (3-1) 

Same procedure is followed to calculate the counterterms of II and 4> form the self energy diagrams. 



3.1 One loop correction of the lepton doublet l L 

The one loop self energy diagram of II consists of the following diagrams |15] 
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Figure 4: One loop self energy diagrams of II- 
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W a 

Ipk 

p + q 

(IV) 



'~:.l 



Taking the divergent part of the above diagrams and following Eq. (13. It) we find the expression for 
the counterterm as 



AZi 



1 



L,l 



16tH 



y v y} + y e y e t + -i B g\ + -i wg \ 



(3.2) 



In the above expression, Y e is the Yukawa coupling coming from the I — e — <f> vertex. g\, g<i are 
U(l)y and SU(2)i gauge coupling constant. £b, £,w are the gauge fixing parameters for the U(l)y 
and SU(2)l gauge fields respectively. Y u Yj term contributes till the heavy field gets integrated out 
from the theory. 



3.2 One loop correction of the Higgs doublet <\> 

The one loop self energy diagram of eft consists of the following diagrams [15 
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Figure 5: One loop self energy diagrams oft 
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Note that the fourth diagram (IV) vanishes in dimensional regularization, hence does not contribute. 
Following the similar procedure as above we find 



AZm 



16vr 2 



2T-~(3-£ B )<7?-^(3-£ w )<?l 



where, in general, 



Tr 



Y U Y} + Y e Yl + 3Y u Yj + 3Y d Y} 



(3.3) 



(3.4) 



Y u , Yd are the Yukawa coupling arising from the q — u — (ft, q — d — <p vertices respectively. Below 
the mass scale M, the heavy field N gets decoupled from the theory. In that case, diagram V 
will not contribute to the self energy of (/>, consequently one needs to omit the term Y v Yj from 
the expression of the T. Here we are not writing down the higgs mass counterterm because it is 
irrelevant for the present work. 



3.3 One loop correction of c (6) up to 0{l/M 2 ) 

Below we give all the one loop diagrams which contribute up to 0(1/ M 2 ) towards the correction 
of c^>. In the following, we proceed with the kind of vertex in Fig. I2N . 



^Taking the other vertex, Fig. [T] will give same /3-function. 
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Figure 6: Diagrams for one loop correction of c^' . 

One loop correction of c^'c*^' which is of the order of O (l/M 2 ) contributes towards the correction 
of c( 6 ' (diagram III). There will be six additional gauge boson diagrams with B replaced by W a . 
Collecting the divergences from all diagrams and with the help of (13. ip . we find the expression for 



Ac i as below 






Acf? = 


1 
16vr 2 


2 (>y e y e t - 




+3 (fo ~ ~) 9^' 



2 (>y e y e t + Y e Y^ Q) ) + \c^c^) + (^ _ 5i y 6) 



In deriving Eq. (J3.5D . we have used the on-shell condition for the external fermions. 



(3.5) 



4 /3- function of c^ 



Now we are ready to compute the /3-function of the parameter, c^' . We collect different expressions 
from the earlier section. From Eq. (|2.9p . the bare non-unitary parameter c B is related to the 
renormalized one as 



.(6) 



To one loop order, Zi L , Z^ and c*- 6 ** from (|3.2p . (|3.3p . (|3.5p read 

1 



i _i 
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A^r„l 



AZ^ 



Ac ( ? 



16vr 2 

1 
16vr 2 

1 
16vr 2 



Yetf + \(LB9l + \iB9\ 

^-\{?-iB)g\-\i?-iB)9l 

2 ( c ( 6 )y e y e t + y e y e t c ( 6 )) + ^"V® + (^ - 5l V 6 ) 

?/.(8) 



+ 3 I iw ~ g I ''2'' 



(4.1) 
(4.2) 

(4.3) 



As we discussed earlier, c^> gets generated below the mass scale M, when the heavy particle gets 
integrated out effectively implying that the term Y u Yj can be dropped from the above expressions. 
Also the expression for T becomes 



Tr 



y e y e t + 3y u yJ + 3y d y c 



(4.4) 



To compute the /3-function we shall use the general equation for /3-function, derived in [14] . The 
/3-function for the non-unitary parameter, c^\ works out to be 



1— loop 



P c (6: 



1 



16vr 2 



- (>)y e y e t + yy e t c ( 6 )) - 2T C ( 6 ) + -c^c*^ + 3 g 2 2 c^ 



(4.5) 



This is the one loop /3-function of the non-unitary parameter which governs the RGE of c^'. It 
is almost similar to the /3-function of dimension-5 parameter, widely known as k in the literature 
[lUCE!]. The term c (5) c* (5) is additional hertH 

5 Summary 

In summary, we have discussed the origin of non-unitarity at high energy for type-I seesaw. Below 
the mass scale of the heavy particle, the dimension-6 Dirac type effective operator leads to the non- 
unitarity of the lepton mixing matrix. We have calculated the one loop correction and derived the 
/3-function which governs RG evolution of the non-unitary operator. Since non-unitarity becomes 



'Coefficient of this term differs from an eariier paper 



significant in TeV scale models, running is low as compared to the operators originated at a higher 
scale. But in fine tuned situation this running may give significant modification to the result. RGE 
of non-unitary operator may also become relevant for precision measurement of neutrino mixings. 

Acknowledgment 

I sincerely thank S. Goswami for inspiring me to do this project and careful reading of the draft. 
I also thank N. Mahajan, S. Maji, S. Mohanty, S. Roy for useful discussions. I am also thankful to 
J. Kersten, A. Raychaudhuri for their feedback. 

6 Appendix A 

Relevant Feynman rules, those are used to compute above diagrams. Arrows on the particle line 
indicate charge flow. Arrows by the side of particle lines indicate momentum flow. Grey arrows 
indicate orientation, introduced in |16j . 
propagators: 

*■ *, , , {B T1- ( f T 1 .-^ ±^t x . 

p 2_ m 2+ iF °Pa / \/\/\/\/\/\/\/ N I -2^- OabOp a 

N [ ~*~ N J i(f+M 3 ) c foci hi « 

p 2_ M 2 +i£ 0ji ' j?+te°ji&Pa 

V p 

when momentum is opposite to the orientation, p should be replaced by —p. f denotes Dirac 
fermions. In case of singlet, the SU(2) indices (i.e. a, (3) should be omitted. 
S 4 interaction: 



d>5 ^ 

\ / 

\ / 

\ / 

\ / . 

/ \ 

• \ 

/ \ 

A / X 

(PR 

gauge boson interaction: 



N AAAAAA/W -h^' 2 9i (p»+p'») s fia N WVWWW -& e/2 9i {p,+p'») rf a 



l -il 




%H e/2 gidp a 6 ji 'y IJ ,PL 
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■&* ,2 920 a pJjil* P L 




%fj, e/2 giS / 3 a 5 j i'Y li PL 



^ e/2 92^J J a fl P R 



0,/ 




^ e glg^8 alj 8 ab 
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o,/ 




2^ 9\9 °<*p 




2^ 51923^ CT a/ = 



Yukawa interaction: 



A'; 



A r , 




-z^/ 2 (r;v. e Q/ ,p L 



' (W 




-i^HY^eaflPL 



N; 
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-W e/2 (Y»)ij tccpPR 




-W e/2 (^)ij CaflPfl 




-W e/2 (Y e %S aP P L 




-i^r< (y*).. 5 aP p L 




-i^/ 2 (y e ) w s aP p R 




-i^ (y e ).. ^p fl 



dimension-5 effective vertex, c' 5 ^: 




c (5) |M £ c*f ) (e 7 5£ Q /3 + e 7 /3£aA-) Pl 




M E c, J (e T 5£a^ + e 7 /3ea,5) Pr 



dimension-6 effective vertex, c^': 
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counterterm of II and 4>: 



hi 



^(AZi^PlSpc 



P P 

c^' counterterm: 
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7 Appendix B 

Here we give the divergent part of different diagrams 
Self energy diagrams of II, Fig. 2} 



(I) 



(II) 



(III) 



(IV) 



16vr 2 



(y v yJ) 5 ia fP L - + UV finite 
\ / ji e 



-^ (y e y e t) .. ^ a fP L X - + UV finite 



i g 2 1 

zrr-^-z-SjiS^BfPL- + UV finite 



i 3 1 

TTT^^glSji^iaZwtfPL- + UV finite 



Self energy diagrams of <p, Fig. [5J 

(I) 



1 - ,2; 1 



„-3Am 5bo — h UV finite 
16n 2 M <Ll 



(II) 



(III) 



(V) 



(VI) 



(VII) 



(VIII) 



i a 2 1 

■Sfia [("3 + Cb)p 2 ~ i B m 2 ] - + UV finite 



167T 2 2 



]s%8p a [(-3 + i w ) p 2 - iwm 2 } i + UV finite 



16vr 2 2* 



^2 (Y v ) jh (yj j ^ ^ [p 2 - 2M 2 ] - + UV finite 



2Tr (Y e YA 5 Pa p 2 - + UV finite 



16vr 



16vr 2 



6Tr (yyj) fyap 2 - + UV finite 



^6Tr (y^) Spap 2 ^ + UV finite 



One loop correction of c^\ Fig. [6l 

(I) 



(II) 



(HI) 



(IV) 



16vr 2 



16vr 2 



if [c^YeYf] .. e 7 /3e Q a (ft - / ) Pz^ + UV finite 



if (yvic^) . e 7/3 e Q(5 (ft - /J P L i + UV finite 



16vr 2 ' 



j i 



// f cWc'W ) I e 7/3 e a< 5 + 7^«5 7 a (ft - ft, J Pl- + UV finite 



1 ,,e9l(6) 



16vr 2 



li e ^cfe lP e aS (ft - /J P L £ B - + UV finite 



Corresponding W boson diagram is three times the P diagram with g\ replaced by 52 • 

(V) 



1 ..e9i(e). 



^ T C S e -yP e <*s (ti ~ /J p ^~ + uv finite 



1 



16tt 2 ^ 4 

Corresponding W boson diagram is three times the B diagram with g\ replaced by 52- 

(VI) 



1 e9f (6) 

1^ -J C ji ^^ 



+ 

12 



P L - + UV finite 
e 



Corresponding W boson diagram is identical with the B diagram with g\ replaced by 52- 
(VII) 



1 ,,e9il&). 



1^ ~l C 3i e 7/3 £ «<5 



1+Cb)^ + ^{ 



P L - + UV finite 
e 



Corresponding W boson diagram is identical with the B diagram with g\ replaced by 52- 
(VIII) 



i ,9i(6), . f . ,\ - 1 



J^2^^ c ji^^B [ft - f/J Pl- + UV finite 

Corresponding W boson diagram is identical with the B diagram with g\ replaced by 52- 

(IX) 

1^2 M e (-f ) 4'V^ (a +4) ^4 + uv finite 

Corresponding W boson diagram is identical with the B diagram with g\ replaced by 52- 

(X) 

~ 16~~2 ^ 8 ji ^ aS P s + € "tP taS ) [Pi + P'i) Pl ~ + UV finite 

In the above expressions, a term proportional to 5~ a 6ss arises in cases of diagram III and all the 
W boson diagrams. We used the identity 

to change 5 ia 5p$ to e 1 p€ ci $. The other part {i.e. dpaS^s) is orthogonal to the non-unitary operator 
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